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In this paper we treat the time evolution of unitary elements in the N level system 
and consider the reduced dynamics from the unitary group U(N) to flag manifolds of 



the second type (in our terminology). Then we derive a set of differential equations of 
matrix Riccati types interacting with one another and present an important problem 
on a nonlinear superposition formula that the Riccati equation satisfies. 

Our result is a natural generalization of the paper Chaturvedi et al (arXiv : 
0706.0964 [quant-ph]). 
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1 Introduction 



In this paper we treat a finite quantum system (N-level system) and consider its unitary 
evolution in detail. 

The unitary evolution for a time-dependent Hamiltonian H(t) (which is not so artificial) 
is given by 

i%^-U(t) = H(t)U(t), U(0) = E N 

where H(t) G H(N, C), U(t) G U(N) and is the identity matrix. In the following we set 
h = 1 and write 

ill(t) = H(t)U(t), U(0) = E N (1) 

for simplicity. 

Now, we consider a reduction of symmetry. Namely, symmetry happens to reduce from 
the unitary group U(N) to a subgroup H. Then the evolution on U(N) "split" into evolutions 
on the subgroup H and corresponding homogeneous space U(N)/H. 

We are interested in the evolution on U(N)/H with special H and call this the reduced 
dynamics from U(N) to U(N)/H for short. 

Especially, for H = U(m) x U(n) (m + n = N) the homogeneous space is the Grassmann 
manifold and the matrix Riccati equation appears naturally, see jT] and [2]. 

Interesting enough, we meet the matrix Riccati equation(s) in several fields in Mathe- 
matics or Physics. In the special case m = n = 1 we have usual (complex-valued) Riccati 
equation. This equation has a mysterious formula called a nonlinear superposition one, see 
the text in detail. 

In the paper we generalize the Grassmann manifold to some Flag manifold. Namely, for 
H = U{1) x U (m) x U{n) (I + m + n = N) we call U (N) / H the flag manifold of the second 
type (temporarily). Then a set of differential equations of matrix Riccati types interacting 
with one another are obtained. This is the main result and we also present a problem on a 
nonlinear superposition when I = m = n = 1. See [3] as a general introduction to this topic. 

As a result it may be concluded that "Riccati structure" appears naturally in the process 
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of symmetry reduction from the unitary group to some subgroup. 



2 Reduced Dynamics on Grassmann Manifolds : Re- 
view 

We set m + n = N (m, n G N) and review some results from pQ, [2] in the case of Grassmann 
mnnifolds. As a whole introduction in this section [3] is recommended. 
The Grassmann manifold is the set of complex vector spaces defined by 

G m , n (C) = {V C C m+n I dimcV = m}. (2) 

Then it is well-known that 

G m>n (C) S U(m + n)/U(m) x U(n) (3) 

and moreover 

U{m + n)/U{m) x U{n) ^ GL{m + n)/B + (4) 
where _B + is the (upper) Borel subgroup of GL{m + n) given by 

B= \l*i * J G GL(m +n) | *i G GL(m; C), * 2 G GL(n; C) 
[ \ * 2 / 

In order to obtain the element of [/ (m+ri) / U (m) x ?7(n) from an element in GL{m+n) / 
it is convenient to use the orthonormalization (method) by Gram-Schmidt. 
For the matrix 

E m 



Z E n 



G= ( I G GL(m + n)/5 + (5) 

where Z G M(n, m; C) we set 

Z I V E n 
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For {Vi, V 2 } the Gramm-Schmidt orthonormalization in the matrix form reads 
V x = V^VlV-i)- 1 ' 2 ==> Pi = ViV} : projection 

V 2 = (Err^-PjV*, V 2 = V 2 (VjV 2 )- 1/2 . 



Explicitly, 



Vi= \ V 2 =( * \ M-^ 

Z \ E n 



where 

L z = E m + Z^Z, M z = E n + ZZK 
Therefore we obtain the unitary matrix 

E -Z 1 " \ ( L~ l/2 
V=(V U V 2 )= I Z I G C/(m + n). (6) 

Z E n I \ M~ 1/2 



Next, if we consider a transformation 









where 


\ u 2 ) 




\ u 2 j 



G — ► G where | £ U(m) x U(n) 

\ U 2 j \ U, 

then the resultant unitary matrix V also transforms like 

Ui 



V — > V 



u, 
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As a result, the procedure to obtain unitary matrices is covariant under the subgroup 
U{m) x U(n), so we can consider V G G mtn . 

A comment is in order. As another definition of the Grassmann manifold the following 
one in terms of projections is also well-known. See [5] and [6] as an elementery introduction 
and [7] as an advanced one. 

{P G M(m + n- C) | P 2 = P, P ] = P, tiP = m} = U(m + n)/U(m) x U(n). 
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The correspondence (V — > P) is as follows : 

-l 



p r I E m \ / — Z< \ / -Em \ / E m —Z^ 

I Z E n \ \ Z E n 



Next, we treat the time-dependent Hamiltonian (which is not so artificial) like 

/ H^t) V\t) \ 

H = H(t)=\ G % + n;C). (7) 

\ V(t) H 2 (t) ) 

Then the (reduced) evolution equation 

i^-V = H(t)V where V = V{t) <^> Z = Z{t) (8) 

reduces to a matrix Riccati equation 

iZ = V + H 2 Z - ZH 1 - ZV ] Z, (9) 

see PP, Note that we ignored the time evolution on U\ and U2, which is out of interest 
at the present time. 

Especially, in the case of m = n = 1 

/ h^t) v(t) \ 
H — \ (10) 

1 I 1 -ztf) \ 

V = . (11) 
y/i + \m 2 \z(t) 1 J 

the (matrix) Riccati equation becomes 

iz — v + {h 2 — hi)z — vz 2 . (12) 

It is very interesting to note that this equation satisfies a mysterious formula called a 
nonlinear superposition. See for example [8] and its references. As a whole introduction 
to the Riccati equation see [9]. 



Namely, let zi, z 2 , z 3 be three different solutions and z be any solution. For the cross- 
ratio defined by 

( , _ z - z 1 z 2 - z 1 

(z, z u z 2 , z 3 ) = (13) 

Z — Z 3 Z 2 — Z 3 

it is easy to see 

z—(z,z 1 ,z 2 ,z 3 )=0 => (z,z 1 ,z 2 ,z 3 ) = k 

with constant k (g C). From this we can express z in terms of k and three solutions z\, z 2 , z 3 
like 

z = kz 3 (z 2 - gi) - z x {z 2 - z 3 ) 
k(z 2 - zi) - (z 2 - z 3 ) 

This is called the nonlinear superposition formula for the Riccati equation. Concerning 
this formula we have the following 

Problem Can this nonlinear superposition formula be derived (or reduced) from usual one 
in the linear equation ? 

3 Reduced Dynamics on Flag Manifolds 

In this section we generalize the results in the preceding section based on Grassmann mani- 
folds to Flag manifolds of the second type (in our terminology). See [TU] and [TT] as a good 
introduction. 

For / + m + n = N (I, m, n G N) the flag manifold of the second type is the sequence of 
complex vector spaces defined by 

*W(C) = {V C W C C l+m+n I dim c V = I, dimcW = I + m}. (15) 

Then it is well-known that 

Fi,m,n(C) = U(l + m + n)/U{l) x U(m) x U(n) (16) 

and moreover 

U(l + m + n)/U{l) x U(m) x U(n) = GL(l + m + n)/B + (17) 



where B + is the Borel subgroup given by 
/*i * * \ 



B 



+ 



G GL{1 + m + n) | *i G GL(Z; C), * 2 G GL(m; C), * 3 G GL(n; C) 



*2 * 

\ o o * 3 / 

Similarly in the preceding section we consider the matrix 



F = 



Ei 

X E m 



\ 



G GHl + m + n)/B+ 



(18) 



\Y Z E n J 



and set 



V 1 



X 
\ Y J 



, v 2 



( ^ 

\ z / 



, v 3 



( o \ 



\E n ) 



For {Vi, V 2 , V3} we perform the Gramm-Schmidt orthogonalization in the matrix form 



like 



Vi = V^VlVi)- 1 ' 2 =>- Pi = fiVi 1 " : projection 

^2 = (^j+m+n - Pi)^ 2 , ^2 = ^(V^) -1 ' 2 =}► P 2 = V 2 V% \ projection 

V 3 = (E l+m+n -P 1 -P 2 )V 3 = {E l+m+n -P 1 )(E l+mHl -P 2 )V 3 , V 3 ^V 3 (V 3 %)- 1/2 . 
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Explicitly, 



V 1 = 



Vo = 



X A~ 1/2 , 

V y J 

( 



(L - (JA-Mt) 



irh-V2 



-A- : yt + a- 1 ^ (L - (JA-Mt) 1 (z - yA-^t)i 



\ 



where 



-XA-!yt - (E m - XA-W) (L - SA-W) 1 (Z - YA^S^ 
K E n - YA- l Y^ -{Z- yA- 1 ^) (L - SA-H^y 1 (Z - YA^S^ / 

(e„ - YA- l Y ] -{Z- YA- l 5 ] ) (L - tfA -1 ^) -1 (Z - yA" 1 ^)^ 
a = Ei + x*x + yfy, L = E rn + Z^Z, 5 = X + Z ] Y. 



x 



-1/2 



Therefore we obtain the unitary matrix 



V = (Vi,V 2 ,V 3 ) e C/(Z + m + n). 



(19) 



Here we list a decomposition of V 



V = 
( 



Ei 

X £ m 



Y Z E, 



a J 



' Ei -A" 1 ^ -A- l Y^ + A- l 5^(L-5A- 1 5^- l (Z -YA' l 5^)^ 
E m -(L- < 5A- 1 (5t)- 1 (Z-FA- 1 <5t)t 
K 



A-i/2 



(L - (5A _1 ^) -1 / 2 



(£7 n - FA-^t _ _ yA- 1 ^) (L - JA" 1 ^) 1 (z - YA^S^)' 1 / 2 
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Compare V with the corresponding unitary matrix in [10J and [1X\ where / = m — n = 1. 



If we consider a transformation 
F — F U 2 



where 



\ 



Uo 



e U(l) x [/(m) x U(n) 



then the resultant unitary matrix V" also transforms like 



V — > V 



u 2 

\ '-J 



As a result, the procedure to obtain unitary matrices is covariant under the subgroup U(l) x 
U(m) x U(n), so we can consider V G F\^ n . 

A comment is in order. As another definition of the flag manifold the following one in 
terms of projections is also well-known. 

{(P, Q) : a pair of projections in M(l + m + n; C) | trP = I, trQ = I + m, PQ = P} 
= U(l + m + n) /U{1) x U(m) x U(n). 



The correspondence (V — ► (P,Q)) is as follows 



/ Ex 



\ 



( F \ 



P = V 







v- 1 = w 







/ 



Ei 



Q = V 



v- 1 = W 
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where 



W 



( Ei -A- 1 ^ -A^yt + A-i^t (L-SA^S^y 1 (Z -YA-^y \ 

X E m -XA- 1 8^ -XA- 1 ^ -(E m -XA- 1 5^(L-5A- 1 5^y 1 (Z-YA- 1 5^y 
\Y Z - YA' l 5^ E n - YA-W - (Z - VA" 1 ^) (L - 5A -1 ^) _1 (Z - FA~ 1 5 t ) t y 



Next, we treat the time-dependent Hamiltonian like 

' ^(t) V^(t) Vj(t) ^ 

H = H(t)= Vl (t) H 2 (t) V*(t) 

\ V 2 (t) V 3 (t) H 3 (t) J 

Then the (reduced) evolution equation 
d 



(20) 



i— V = H(t)V where V = V(t) 

(Jjv 



> X = X(t), Y = Y{t), Z = Z(t) 
gives a set of matrix Riccati equations interecting with one another 



(21) 



iX 
iY 
iZ 



V x + H 2 X - XE X - XV}X + V3 y - XV^Y, 
V 2 + H 3 Y - YHi - YV 2 Y + V 3 X - YVfX, 

v 3 + h 3 z - zh 2 - zvlz + (ZX - Y){yl + v\z). 



This is our main result in the paper. 

Especially, in the case of I — m — n — 1 

( h(t) v^t) v 2 (t) \ 



H 



vi(t) h 2 (t) v 3 (t) 
V v 2 (t) v 3 (t) h 3 (t) ) 

I 1 -(x+yz) xz 



V 





VA1A2 


v 7 ^ 


X 


l-[xz-y)y 


— z 




VA1A2 




y 


z+x(xz— y) 


1 




VA1A2 





(22) 



(23) 



(24) 



where 



Ai = 1 + \x\ 



\y\\ 



I I 2 1 I \2 

\z\ + \xz — y\ 
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the interacting Riccati equations become 

ix = vx + (h 2 — hi)x — vix 2 + v 3 y — v 2 xy, 

iy = v 2 + (h 3 - hi)y - v 2 y 2 + v 3 x - Vixy, (25) 
iz = v 3 + (h 3 - h 2 )z - v 3 z 2 + (xz - y)(v 1 + v 2 z). 

Concerning these equations, our interest is the following 

Problem In this system is there a nonlinear superposition formula like flHj) ? 

This is an important problem on nonlinear superposition. However, we cannot find such 
a formula at the present time, so we leave it to readers as a challenging problem. 



4 Discussion 

In this paper we considered the reduced dynamics from the unitary group U(N) to flag 
manifolds of the second type and derived a set of differential equations of matrix Riccati 
types interacting with one another. 

We also presented an important problem in the special case (/ = m = n = 1) on 
generalization of the nonlinear superposition formula that the Riccati equation satisfies. 

In last, let us make a comment. We can generalize our construction to flag manifolds of 
the third type (in our terminology) 

F k ,i, m ,n = U(k + l + m + n)/U(k) x U{1) x U(m) x U(n) 

where 



F k , lim , n {C) = {WCVCWCC 



k+l+m+n 



dim c W = k, dim c V = k + l, dim c W = k + l + m}. 



In fact, starting from the matrix 
/ 



E k 
K E x 



L X En 







MY Z E r . 



E GL(k + l + m + n)/B^ 
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we can trace the same process as the text. Explicit calculations done are of course very hard, 
see [T2] . 

When k = l = m = n = l some explicit calculations have been done by Picken |10j . 
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